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Code division multiple access (CDMA) is a core technology of today's wireless communication system 
that employs data transmission between multiple terminals and a single base station. In the general 
scenario of a CDMA system, the binary signals of multiple users are modulated by spreading codes 
assigned to each user, and these modulated sequences are transmitted to a base station. The base station 
receives real-valued signals, which are a superposition of the modulated sequences and possible noise. 
A detector at the base station then extracts the original binary signals from the received signals using 
& • knowledge of the spreading codes of the users. This task is termed user detection. 

The Bayesian framework indicates that the optimal detection performance can be achieved by mul- 
tiuser detection, which simultaneously infers the signals of multiple users. Although performing such 
optimal detection exactly is computationally difficult, this difficulty can be approximately resolved by 
employing belief propagation (BP), and the macroscopic behavior of the BP-based algorithm is well 
captured by the natural iteration of the saddle point equations of replica symmetric (RS) analysis [Q. 
Unfortunately, the performance of the BP-based algorithm deteriorates as the discrepancy between the 
true and assumed channel models becomes larger. In particular, this algorithm exhibits dynamic instabil- 
ity when the intensity of the channel noise is set below a certain critical value, which is in accordance with 
the de Almeida- Thouless (AT) instability in equilibrium, analysis J] . This correspondence of two types of 
instability has motivated the development of a multiuser detection algorithm based on survey propagation 
(SP), which is a recently developed algorithm that is inspired by the one-step replica symmetry breaking 
(1RSB) solution of equilibrium analysis |2|. This is the purpose of the present study 2 . 

We focus on a CDMA system that uses binary phase shift keying (BPSK) symbols and K random 
binary spreading codes of the spreading factor N with unit energy over an additive white Gaussian noise 
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(AWGN) channel. For simplicity, we assume that the signal powers are completely controlled to unit 
energy, but the extension to the case of distributed power is straightforward. In addition, we assume that 
chip timing and symbol timing are perfectly synchronized among users. Under these assumptions 0], the 
received real- valued signal can be expressed as 
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where /i G {1, 2, . . . , N} and k S {1, 2, . . . , K} are indices of samples and users, respectively, s^k £ {~ U 1} 
is the spreading code with unit energy independently generated from the identical unbiased distribution 
P(Sfj. k = +1) = P(s /lfc = -1) = 1/2, b° k is the bit signal of user k, is a Gaussian white noise sample 
with zero mean and unit variance, and o~o is the standard deviation of AWGN. 

In the Bayesian framework, the multiuser detection problem can be formulated as a problem to 



infer K bit user-signals b° — (bi) from the posterior distribution P(b\y) 



with 



knowledge of the spreading codes of the users {s^fc}, where P(y fi \b) — -^== exp — ^ (y^ — u M ) 2 
b = (6^) = {+1, —1} K , = -j= Yl,k=i s ^kbk and o is the noise parameter assumed at the base station. 
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2 A similar but different approach to an improved algorithm can be found in 
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For such an inference problem, the replicated Bethe free energy formalism [5] expresses SP as 
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where t represents the number of updates and auxiliary field distributions TT t k ^ ll {h) and 7r„_,fc(/&) are 
termed surveys, the stationary solutions of which characterize the replica symmetric stationary points 



of the replicated Bethe free energy. h({hjjtk}) — tanh / h t — \ anc ^ x cor " 

responds to the replica symmetry breaking parameter of the 1RSB analysis. At each update, the 
marginal distribution is approximately assessed as P(b k \y) — 2& P{b\y) ~ J dhp k {h) 2c osh(h k ) ' wnere 
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Unfortunately, executing this algorithm exactly is practically difficult because the computational 
cost for assessing h({hj^ k }) grows exponentially with respect to the number of users K. However, 
this difficulty can be overcome using Gaussian approximation when K and N are large while keeping 
(5 = K/N ~ 0(1), which leads to a computationally tractable algorithm 
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MI = « t ; . Here, we set Dz = -%==e 2 Z . This provides the maximizer 
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of posterior marginal (MPM) estimate, which minimizes the bit-wise detection error probability in the 
ideal case, as bk — sgn(m k ) at each update, where sgn(u) = u/\u\ for m ^ 0. The computational cost 
for carrying out this algorithm scales as O(NK) per update when the one-dimensional integrations for 
assessing rrik and Mfc are approximately evaluated with a fixed amount of time. 

A tree approximation, which ignores time correlations of variables between successive updates, implies 
that the macroscopic behavior of this algorithm is described by the natural iteration of the saddle point 
equation of the 1RSB analysis. On the other hand, numerical experiments indicate that the discrepancy 
between the trajectories of macroscopic variables and those of the 1RSB saddle point equation is relatively 
large when the AT stability of the RS solution is lost, in particular, during the transient stage. However, 
the macroscopic property of the fixed point excellently accords with that predicted by the 1RSB solution. 
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